Simple Physical Picture of the Overhauser Screened Electron-Electron Interaction 
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As shown by Overhauser and others, the pair-distribution function g(r) of a many-electron system 
may be found by solving a two-electron scattering problem with an effective screened electron- 
electron repulsion V(r). We propose a simple physical picture in which this screened repulsion 
is the "dressed-dressed" interaction between two neutral objects, each an electron surrounded by 
its full-coupling exchange-correlation hole. For the effective interaction between two electrons of 
antiparallel spin in a high-density uniform electron gas of arbitrary spin polarization, we confirm 
that this picture is qualitatively correct. In contrast, the "bare-dressed" interaction is too repulsive, 
and does not have the expected symmetry Vj|(r) = Vj.f(r). The simple original Overhauser model 
interaction, independent of the relative spin polarization £, does not capture the ^-dependence of 
the correlation contribution to g(r — 0). 
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The quantum mechanical many-electron problem is no- 
toriously hard if all its degrees of freedom are taken 
into account. For both practical computational and 
conceptual purposes, however, it can often be replaced 
by a one- or two-electron problem with an effective 
external potential or electron-electron interaction, re- 
spectively. The effective potential that shapes the or- 
bitals of the one-electron problem in Kohn-Sham density 
functional theoryiS has been intensively explored, but 
the effective screened interaction that shapes the gem- 
inals of the two-electron problemi^SiSiLSi has received 
less attention. Here we propose and provide some sup- 
port for a physically-appealing "dressed-dressed" pic- 
ture based upon the interaction between two neutral ob- 
jects, each being an electron dressed by its surrounding 
exchange-correlation hole. In this picture, the "bare- 
bare" Coulomb repulsion 1/r is strongly screened out 
over the Wigner-Seitz radius r s . 

Overhauser 7 showed that the singlet geminals of an 
effective two-electron scattering problem can be used to 
estimate the on-top pair-distribution function g(0) in a 
spin-unpolarized (£ = 0) three-dimensional electron gas 
of uniform density 



(1) 



(We use Hartree atomic units where H = m = e — 1.) 
Overhauser used an effective "bare-dressed" interaction 
between a bare electron and a neutral object composed 
of another electron and a concentric sphere of positive 
background charge of density n and radius r s . Gori- 
Giorgi and Perdew<2ii£ used the same effective interac- 
tion, but solved the Overhauser model exactly and found 
a pair-distribution function g(r) in close agreement with 
that of Quantum Monte Carlo calculations over the whole 
short-range region r < r s , for the physical density regime 
1 < r s < 10. For the high-density (r s — * 0) limit, they 
found good agreement with the exacfeii*^ g(r) to order 
r s . Since then, there have been many related studies 
for the three- or two-dimensional electron gasji^iiiSiiSii 7 . 



sometimes using constructions of self-consistent effective 
interactions following the general "bare-dressed" picture 
of Overhauser Sum rules for the scattering phase shifts 
have also been derived^ 

In this work, we consider a three-dimensional uniform 
electron gas with relative spin polarization 

C=(n T -n i )/n, (2) 

where n = nj + nx i s the total density of Eq. (JJJ . The 
pair-distribution function is then 

g(r) = (^)%ttM + (¥)%UW + ^STiW, 

(3) 

where only g^ contributes at r = because of the Pauli 
principle. ng(r) is the average density of electrons at r 
when an electron is at the origin, and n [g(r) — 1] is the 
density of the exchange-correlation hole at full coupling 
strength, which carries a charge equal and opposite to 
that of the electron it surrounds: 10 



dr 47rr 2 n [g{r) — 1] 



(4) 



with the same equation for n-f [<7jj (r) — 1] and [,9li( r ) ~ 
1]. We focus on the effective interaction T / tl( r ) between 
two electrons of opposite spin in the high-density (r s — > 
0) limit, since in this case correlation can be neglected 
and the interaction is purely electrostatic. Thus we can 
evaluate the "bare-dressed" and "dressed-dressed" mod- 
els exactly and compare the predictions of both to the ex- 
act pair-distribution functio n 11 ! 12 : 19 whose short-ranged 
part is dominated by Vj-j. (r). We do not explicitly discuss 
the electron-electron scattering effects on transport prop- 
erties, which are a second important application of the 
effective two-electron problemAiiSii 7 . We note however 
that the expected symmetry V^i = Vj.j of the effective 
interaction for £ ^ is only achieved by the "dressed- 
dressed" picture, not by the "bare-dressed" one. 
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Spin-unpolarized gas- In the Overhauser approaches 
to electronic correlation in the unpolarized (£ = 0) uni- 
form gas, the many-electron problem is reduced to a 
scattering event between two electrons in a suitable ef- 
fective potential V(r,r s ), with a corresponding radial 
Schrodinger equation: 



dr 2 



£(1 + 1) 



V(r,r 8 ) + k 2 u e = Q, 
u e = krRi(r, k, r s ). 



(5) 



The presence of the other electrons is taken into account 
in two ways: (i) via V(r,r s ), (ii) via an average over the 
possible relative momentum k = ^|ki — k^l of the scat- 
tering event. The exchange symmetry between the two 
electrons is ensured via a proper summation over the par- 
tial waves £] the resulting spin-resolved pair-distribution 
functions are then 9 
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FIG. 1: £ dependence of the high-density (r 3 — > 0) correction 
to the on-top value A(£)/A(0) [see Eq. 1121 1, The result from 
the "dressed-dressed" potential of Eq. 1151 is compared with 
the exact calculation^-*— with the result obtained from the 
original Overhauser potential of Eq. (|5J , and with the scaling 
relation proposed in Ref. 9 (GP). 



9n(r,r s ) = (£(2*+l)i#(r,fc,r,)\ ( 6 ) 
U=o / 

g n (r,r s ) = 2 ( jr (21 + l)R 2 e (r, k,r.)) , (7) 
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where the symbol (• ■ • ) denotes the average over the 
probability p(k) (obtained from the momentum distri- 
bution of the ideal Fermi gas 9 ). Overhauser's origi- 
nal choice 7 for V(r,r s ) was the potential of an electron 
surrounded by a Wigner-Seitz sphere of uniformly dis- 
tributed positive charge: 

V(r,r s )= i(l + 4-§) ( r <r s ) 

(r > r s ), (8) 



where 



r/r s 



(9) 



is a scaled variable. As said, this simple potential gave 
surprisingly accurate results^ for the short-range (r < r s ) 
part of the unpolarized-gas g(r), at metallic and lower 
electron densities. The result for the high-density (r s — > 
0) limit was also quite accurate: the form of the screened 
Overhauser potential ensures that the correction to the 
noninteracting gas for r s — > is of first order in r s , as in 
the exact perturbative resulttii 

<W (s, r s ^0)= g<® ( S ) + r s (s) + o(r s ), (10) 

where g( ' is the pair-distribution function of the nonin- 
teracting gas. (Eq. (|10[1 is valid for r <C \/rs-) In par- 
ticular, for the value of the tl pair-correlation function 
at contact (r = 0), the solution of the Overhauser model 
gives** g-\i(r — 0, r s — > 0) = 1 — 0.694 r s + o(r s ), in reason- 
able agreement with the exact resun— 1 — 0.732 r s +o(r s ). 
Nagy et alM> have shown that the high-density form 
5ti( r = 0, r s — > 0) = 1 — A r s + o(r s ) is guaranteed when 



Eqs. ©-0 employ a screened potential with screening 
length cx r s . For finite r, the r s — > form of Eq. I|10|l is 
satisfied, within the Overhauser approach, if the poten- 
tial V(r,r s ) is such that 



V(r,r s ^Q) = -U{s 



(11) 

fulfills Eq. (TTJ at 



The Overhauser potential of Eq. 
all r s . 

Spin-polarized gas - In the original formulation of the 
Overhauser modelji information on the spin polarization 
state of the electron gas only enters through the probabil- 
ity distribution for the relative momentum k. The poten- 
tial, purely based on classical electrostatic arguments, is 
independent of £■ The probability functions p" a (k) are 
given m Eqs. (42)-(44) of Ref. ©, where, however, the 
calculations for the Overhauser model with ( ^ have 
not been carried out. Instead, a scaling relation has been 
proposed. 

Here, we carry out the calculations for the high-density 
limit with the correct p^ (k) , and we find a very weak £ 
dependence of the first-order correction A(£) to the on- 
top value, 



g n {r = 0, r s -> 0, C) = 1 - A(C) r s + o(r s ), 



(12) 



as shown in Fig.^ This is due to the weak£ k dependence 
of the short-range part of the s-wave radial wavefunc- 
tion Ro(r — > 0, k, r s ) of Eq. (J5J. An explicit dependence 
on £ in the effective potential is thus needed in order 
to reproduce the correct behavio»ii*i*i of the short-range 
part of g(r) in the spin-polarized electron gas. Moreno 
and Marinescui^ have recently applied the Overhauser 
model to the two-dimensional electron gas, finding an 
extremely weak £ dependence of the on-top value. Our 
Fig. ^ suggests that their result could be an artifact of 
their (^-independent effective interaction. 

Effective interaction for opposite-spin electrons - In 
the high-density limit, a simple physically-motivated ef- 
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fective potential for antiparallel-spin interactions, which 
depends on £ and has the symmetry Tl=lt> can be ob- 
tained in the following way. Consider two electrons of op- 
posite spin in a uniform electron gas in the high-density 
limit. Each electron induces around itself an exchange 
hole, forming a neutral object. The effective potential 
can be approximated with the electrostatic interaction 
between two neutral or "dressed" objects. When £ — 0, 
each electron is surrounded by a compact exchange hole, 
leading to effective screening of the Coulomb repulsion. 
But as £ approaches 1, the exchange hole around the 
minority spin will become shallow and broad, so the 
Coulomb repulsion will be less well screened. 
The two charge distributions are then 



Pi( x ) 

P2(x) 



<S(x) 



^ T [ 5 ! T (x)-l] 



(13) 
(14) 



rj(x-r) +n 4 .[^ i (x-r) -1], 
and the corresponding electrostatic potential is given by 



V(r,r s ,C) 



dx / dx 



, pi(x)p 2 (x') 



(15) 



V(r,r s ,Q can be computed analytically: its Fourier 
transform V(k, r s , £) is equal to 



V{k,r s ,C) 
with 

vi = 

V2 = 
V3 = 



4tt 
k 



2 +v 1 (k,r s ,C) +v 2 {k,r s ,C) + v 3 (k,r s ,(), 



[S^(fc,r„C)-l]|f 

4-7T 



(16) 
(17) 



47T 



[SP(k,r a ,0-l][SjiHk,r ai Q-l]jp, (18) 




FIG. 2: dependence of the dimensionless "dressed-dressed" 
potential U(sX) calculated from Eq. 1151 . 




FIG. 3: High-density (r s — > 0) Tl correlation holes computed 
from the "dressed-dressed" potential of Eq. 1151 for different 
values of the spin polarization (. 



behavior is very similar to the exact one of Fig. 1 of Ras- 
solov et (i/pLi This is more evident in our Fig. ^ where 
the function A(£)/A(0) is compared with the exact result. 



where S™ are the exchange-only static structure factors, While the ^ de P endence of »u ^ obtained from the 



C(T(T 3 fc _1_ / k 

D x ~ 4 k" F 16 \k F 

1 



(fc < 2k F ) 
(fc > 2k F ), 



(19) 



with k a F 



(l+sgn^O^fcF, k F = (f) 1 ^ 1 , and 
sgn(cr) = +1 for spin-| and —1 for spin- J. electrons. The 
exchange-only pair-distribution function g x only depends 
on r s through the scaled variable s — r/r s . This ensures 
that V(r, r s ,() = ±U(s, (), as required by Eq. JTTJl. The 
dimensionless potential U(s, £) is screened for s > 1, and 
goes to zero, when s — > oo, as s~ 4 . Its £ dependence is 
the one expected from the qualitative arguments given 
above, as shown in Fig. when £ — > 1 the potential is 
less and less screened; for ( exactly equal to 1 (but only 
in this case) U(s — > oo, ( — 1) goes to zero as s~ 2 . 

Using the effective potential U(s,() in the Over- 
hauser scheme, we calculated the fj, high-density pair- 
correlation functions, 9$, for different values of the spin- 
polarization £. They are shown in Fig. El the qualitative 
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FIG. 4: High-density (r s — > 0) Tl pair-correlation function for 
the C = gas obtained from different screened interactions: 
the "dressed-dressed" potential of Eq. I|15|l . the original Over- 
hauser potential of Eq. (|SJ , and the "bare-dressed" potential 
of Eq. 12UII . The exact calculation of Rassolov et al^ is also 
reported. 
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s = r/r s 

FIG. 5: Comparison of the bare Coulomb potential with dif- 
ferent simple screened potentials for the £ = gas in the 
high-density (r s — > 0) limit: the "bare-dressed" potential of 
Eq. 1201k the "dressed-dressed" potential of Eq. 11511 . and the 
original Overhauser potential of Eq. JHJ . All curves have been 
multiplied by r a . 

simple potential U(s, £) is rather good, the quantitative 
agreement with the exact result when £ = is less ac- 
curate than the result obtained with the original Over- 
hauser potential. This is shown in Fig. 0] we see that 
for small s, obtained with U(s, £) of Eq. (|15|) is too 
deep, while the original Overhauser potential of Eq. © 
gives a result which is slightly less deep than the exact 
one. This means that the original Overhauser potential of 
Eq. JSJ is slightly too screened in the r s — > limit, while 
U(s, C = 0) of Eq. i|15|) is not screened enough in the same 
limit. The "exact" effective potential for the high-density 
limit should thus lie in between the two curves "Over- 
hauser" and "dressed-dressed" of Fig. \5\ In the same 
figure we also show the bare Coulomb potential, and the 
"bare-dressed" potential (obtained from the interaction 
of a "bare" electron with a "dressed" electron, i.e., sur- 
rounded by its exchange hole), whose Fourier transform 
Vi(fc,r s ,C) is 

4-7T 

Vi(k,r s ,Q= ^ +vi(k,r s ,0, (20) 



where fi is given in Eq. (|16(l . The "bare-dressed" po- 
tential is "philosophically" closer to the original picture 
of Overhause»2il and to the high-density limit of the 
self-consistent Hartree approximation of Davoudi et alM> 
We see that the "bare-dressed" potential is much less 
screened that the "dressed-dressed" one and thus corre- 
sponds to a a deeper (i.e., further from the exact result) 

gi 1 ^, as shown in Fig.^1 

The "bare-dressed" potential encounters severe prob- 
lems for the calculation of A(C)/A(0) of Fig. ^ When 
C — * 1, each majority | electron dresses itself in an ex- 
change hole deeper and more short-ranged than for £ = 0, 
while each minority j electron undresses. So the interac- 
tion between a bare J, and a dressed f becomes less re- 
pulsive as £ increases from 0, reducing A(£)/A(0). If we 
try to symmetrize using the interaction of a hypotheti- 
cal bare | with a dressed J,, we find that this interaction 
tends to the unscreened 1/r as £ — * 1. 

Conclusions - We have proposed a simple "dressed- 
dressed" picture for the effective screened electron- 
electron interaction that shapes the geminals and thus 
the pair distribution function of a many-electron sys- 
tem. In this picture, the interaction is between two neu- 
tral objects, each an electron dressed by its exchange- 
correlation hole. For two electrons of opposite spin in a 
high-density electron gas of arbitrary spin polarization, 
where the "dressed-dressed" and "bare-dressed" interac- 
tions can be evaluated exactly, we have shown that the 
"dressed-dressed" picture is qualitatively correct. In fu- 
ture work, it may be possible to construct the "dressed- 
dressed" V aa ' (r) for all r s and £, using density functional 
theoryi£ to describe the additional exchange-correlation 
terms that arise when a' — a or r s 0. 
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